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, A systematic analysis of the complex lineshape at (7r,0) of heavily overdoped Bi2Sr2CaCu208+5 

' is presented. We show that a coherent component in the quasi-particle excitation is correlated with 

■ the superfluid density throughout the entire doping range and suggest it is a direct measure of the 
CN . order parameter of high temperature superconductors. 
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The off-diagonal nature of the order parameter of a superconductor makes it difficult to be 
measured directlyf])]. For conventional superconductors, the appearance of an energy gap in 
the single particle excitation spectrum characterizes the superconducting state; as a result, 
the observation of the gap is virtually equivalent to the observation of the superconducting 
order parameter. On the other hand, for high temperature superconductors (HTSC's), 
because of the pseudogap phenomenon in the normal state of underdoped systems, such a 
strict link between the gap and off-diagonal long range order (ODLRO) breaks. There have 
been suggestions between ODLRO and the (vr, vr) resonance peak in the dynamic structure 
factor measured by inelastic neutron scattering experiments around 40meV in HTSC's P|. 
However, because a significant amount of spin fiuctuations exist in this region and evolve with 
temperature in the normal state of underdoped systems, a direct and clean measurement of 
the superconducting order parameter is still not achieved. 

We pointed out earlier that in the single particle excitation spectrum a sharp peak, 
which develops near (tt, 0) of the Brillouin zone below the superconducting transition 
temperature (Tc), exhibits similar behaviors as the superfiuid density j^. This so-called 
superconducting peak (SCP) is observed by angle-resolved photoemission spectroscopy 
(ARPES) in many HTSC's, first on Bi2Sr2CaCu208+5 (Bi2212), later on FEaaCugOT.^, 
and Bi2Sr2Ca2Cu30iQj^s[4:\- Being a probe for single particle excitation, ARPES cannot 
measure collective properties such as the superfiuid density (or phase stiffness) directly. 
However, it could be measuring condensate fraction, a single particle quantity closely related 
to superfiuid density. Condensate fraction is defined as fraction of particles participating 
in pairing, and thus is the direct measurement of the ODLRO. Unlike the gap, the correla- 
tion between the SCP and superconducting order parameter works well in the underdoped 
regime. However, its validity in the overdoped regime has been questioned. Some claim that 
the SCP increases upon overdoping, while the superfiuid density decreases j^. Others claim 
that there is no SCP on the overdoped sample, and the peak-dip-hump (PDH) lineshape 
associated with the SCP is simply due to the bilayer band splitting (BBS) effect 0, Q, 0, 0] ■ 

In this paper, we present a global study on the photon energy (hu) and temperature 
(T) dependent ARPES spectra of heavily overdoped Bi2212 at (vr, 0). We show that there 
is a substantial SCP in the heavily overdoped regime and its correlation with the ODLRO 
holds across the entire doping range. This also gives a unified understanding of the PDH 
lineshape in the BBS circumstancejlO]: the anti-bonding band (AB) and bonding band (BB) 
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each develop their own PDH hneshapes by transferring parts of their spectral weight into 
the newly-created coherent quasiparticles (AB peak and BB peak, respectively). 

ARPES experiments were performed at beamline V-4 of Stanford Synchrotron Radiation 
Laboratory (SSRL). Data were collected with He-I light from a He discharge lamp (for the 
T-dependence) and polarized synchrotron light from a normal incidence monochromator (for 
the /ii/-dependence) , where the contamination to the spectra from the second order light is 
extremely weak, ideal for a quantitative comparison between spectra taken with different 
/ii/'s. Detailed descriptions of experimental settings and sample preparation can be found 
in Ref.Q. 

The existence of the coherent peak has been discussed extensively for optimally doped 
and underdoped systems based on doping and T dependence Here we recap this with 
the data taken on an optimally doped Bi2212 (Tc =90K) sample in Fig. 1. As shown in Fig. 
la, the sharp coherent peak grows out of a smooth spectrum with decreasing temperature 
at (vr, 0). Fig. lb shows how this peak behaves along F — M. We found that the SCP 
is strongest when the hump is within 120~150 meV below Ep'-, away from this region, the 
SCP rapidly dies out while dispersing slowly with the hump to higher binding energies. The 
extracted superconducting state dispersions of the peak and hump in the main band are 
plotted in Fig. Ic. We emphasize that the SCP is not the AB, otherwise, one would obtain 
a bilayer splitting increasing towards (0, 0) (Fig. Id), while band theory as well as symmetry 
arguments gives zero splitting along (0,0) — in. 7r)|ll|. 

The complication in the heavily overdoped regime originates from the resolution of the 
AB and BB in the normal state, which forms a PDH-like lineshape with the sharp AB 
in the lowest binding energy position[10]. Kordyuk et al. have implemented a model fit 
of two bilayer-split bands to study the /iz/-dependent superconducting state (vr, 0) spectra 
in this doping regime|3]. A qualitative agreement between the fitted /iz/-dependent matrix 
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elements of the peak and hump and those of bilayer-split bands predicted by theory 
was argued to justify the peak and hump being due to the AB and BB, at least in the 
overdoped regime^]. However, we found that although the fitting function is able to fit 
each single spectrum well at different hu^s, it faces difficulty when more physical constraints 
are included. One example is shown in Fig. 2a. These spectra are taken at (vr, 0) on 
an overdoped =65K sample (OD65) at different temperatures. We start with the nine 
normal state EDC's to obtain the global parameters through a global fit (described below). 
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The fit gives an excellent agreement with the experiments which confirms the qualification 
of the chosen form of spectral function for describing the normal state features and, notably, 
the thermal sharpening-up process on the lineshapes upon T decrease. However, with these 
parameters, the discrepancy between the eight simulated and experimental superconducting 
state spectra increases upon going further into the superconducting state. This demonstrates 
a drastic lineshape variation beyond the description where only two bilayer-split bands are 
considered. The much more rapid intensity accumulation on the experimental peak feature 
below Tc requires further investigation. 

Reflecting the behavior of the SCP in underdoped, optimallv doped, and even shghtly 
overdoped samplesQ], we extend the fitting function in Ref.js] by adding two spectral 
components to account for the rapid spectral weight accumulation near Ep. 

I{Lj,T,hu) = Io{T,hu) ■[{J2Ja{uJ,T,hu) ■ f{uj,T)) 

a 

miuj, T'ihiy)) + B{iu, T)] + h{T, hu) 

, where / is the Fermi function, R the /iz/-dependent resolution Gaussian, B the momentum- 
independent but T-dependent empirical background function obtained by a separate fit on 
the featureless energy distribution curves (EDO's) around (vr/2, vr/2)[l^, Jq and Ji the linear 
intensity coefficients. The summation is over the spectral intensity of AB hump, BB hump, 
AB peak and BB peak, given by, respectively, 

Jah{uJ, T, hu) = Mah{hv) ■ Ca(T) ■ Ah{uJ, T, a, Eah), 

Jbhi^, T, hv) = Mbh{hv) ■ Cb{T) ■ Ahiuj, T, a, Ehh), 
Japiu;, T, H = Mapihu) ■ (1 - Ca{T)) ■ Ap{u;, Ta{T),eap), 
Moo, T, hu) = Mbp{hu) ■ (1 - a(T)) ■ Ap{uj, Tb{T),ebp) 

, where M is the matrix element, C the remaining spectral weight in the hump, A^ and Ap 
the self-normalized spectral functions for the hump and the peak, respectively, given by 





We adapt the empirical J2"{^,T) = J {auy + {j3TY from Ref.Q, where a and j3 are made 



globally fittable {hv- and T-independent) to reflect the doping and impurity level of samples, 



and ^ is a normalization factor. T{T) is a T-dependent linewidth (FWHM). Note that three 
major extensions in the fitting function have been made. First, two spectral components 
for the SCP's are added w.:.Ke. Gua.,a„ Un.Kape. iu.,fled .^e STM ....Q; 
second, spectral weight transfer (SWT) (= 1 — Ca(b),0 < Ca(b)(T) < 1) is proposed for the 
formation of SCP, which is a pronounced character of many correlated systems and is also 
based on the empirical observation of the fact that the integrated spectral weight is roughly 
T-independent in the (tt, 0) region ^^1; third, a realistic SCG opening extracted with the 
leading edge gap method on the same sample introduces an effective T-dependent energy 
shift to the bilayer-split bands in the superconducting state. 

In order to retrieve reliable quantitative information from a multi-dimensional fit to a 
set of interconnected EDO's, we developed a non-linear least-squares fitting routine, which 
has been applied in the normal state fit in Fig. 2a. The goal is to achieve a systematic 
understanding of both the hv- and the T-dependent behavior of the lineshape. Physical 
constraints are strictly followed here, for example, in the the /iz/-dependent set, all functions 
of T are global variants to all EDO's at the same T while the functions of hv (matrix 
elements) are locally specified, and in the T-dependent set, vice versa. The renormalized 
band energy position e is always a global fitting parameter for the EDO's of samples at 
the same doping level. Without any manual tuning in the global parameters, the new 
fitting function can yield a robust, physically-constrained global fit for a large set of EDO's 
regardless of the initial values input. 

It has been shown before that the intensity of the BB peak drops quickly when the BB 
disperses deep below the Ep at (tt, 0) jl^. For simplicity, we first neglect Jhp term, and 
achieve an excellent global fit for all the spectra shown in Fig. 2b-d simultaneously, which 
include the aforementioned T-dependence data and the /izz-dependence (vr, 0) data of a Pb- 
doped OD65 (Pb-OD65) taken above and below T^.. In Fig. 3a, the T-dependent energy 
positions of the maxima of three spectral features are shown. In the superconducting state, 
the AB peak is located quite close to the AB hump, which makes it very difficult to resolve 
them separately in experiments. For the AB (BB) hump maximum, the slight decrease 
(increase) above T^ is fully determined by the T-dependence of the self-energy with the 
fitted a = 0.715 ± 0.006eV^~^ and P = 2.57 ± 0.03eV~^ while the superconducting state is 
mainly affected by the opening of the gap. The spectral maximum shift in the BB hump 
is T-insensitive as compared with the AB hump; for the AB peak which turns up in the 
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superconducting state, its almost-constant location (—18.01 ± O.OSmeV) seems to define 
the maximum SCG opening at (vr,0). This contradicts the scenario where the peak near 
Ep is just the AB which should change its position with the SCG opening. The extracted 
T-averag ed BBS amplitude is 84.39 ± 0.07 mevQ, 

consistent with Ref.llOH. 
The fitted matrix elements as a function of hu for the three spectral features are shown 
in Fig. 3b. Matrix elements of AB and BB hump vary in the same range of magnitude but 
with an anti-phase-like /ii/-dependence as predicted by theoryjl^, Q- Note that the self- 
normalization of the spectral function is required not only for a well-defined SWT analysis 
but also for a physically comparable matrix element between different spectral components. 
< Mf,h > / < Mah >= 0.966 ±0.007, while the matrix elements of AB peak shows a different 
/ii^-dependence with those of the humps and (< Map > / < Mah >= 0.610 ± 0.007). This 
deviation is not unexpected because SCP is of redistributed spectral weight related to the 
superfluid, whose wavefunction may have different properties with the hump counterpart. 

The T-dependent SWT from the above fit is summarized in Fig. 4a. The most prominent 
feature is the superconducting peak emerges only below Tc, while the SCP at lower doping 

. This may associate with the super- 



2l|. The SWT of AB is extrapolated 



could emerge lOK above Tc (see Fig. la and Ref. 
conducting fluctuations in the lower doping regime 
to be 0.223 ± 0.028 at K, a value far below 1, which casts doubt on the scenario where the 
entire lower- lying AB peak possesses almost all the weight of AB quasiparticle 2M- 



present the result of a 4-component fit with both Jap and Jbp terms. The BB peak weight 
extrapolates at K to be only 0.015 ± 0.008, much smaller than that of AB, which also 
justifies the 3-component fit. 

In Fig. 4a, the superfluid density measured by muon spin relaxation (/iSR) at the same 
doping level is also displayed, which resembles the SWT in the T-dependence. In Ref. jsj, 
we proposed a phenomenological substraction of a smooth background to retrieve the SCP, 
which was further quantified by defining the so-called superconducting peak ratio (SPR) as 
the SCP intensity divided by the overall spectral weight in (vr, 0) spectrum taken at hv=22.1 
eV and low temperature. Although the SWT is more physically well-defined than SPR in 
that it is free of the matrix element effects and takes into account the multi-component 
lineshape, when the hump is broad near or below optimal doping, SPR obtained by the 
phenomenological procedure reasonably characterizes the SCP. Thus, we calculated SPR for 
OD65 at hu=22.7 eV at lOK based on the globally fitted AB peak. It is plotted in Fig. 
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4b together with SPR's at lower doping levels reproduced from Ref. P]. One can clearly see 
that SPR decreases with doping in the heavily overdoped regime, similar to the superfiuid 
density [3, 23 1. To account for this decrease, it has been suggested that the overdoped system 
is phase-separated into superconducting regions and hole rich normal metal regions [23]. In 
any case, these two intimate correlations between the SCP and superfiuid suggest the SWT 
is a direct measurement of the superconducting order parameter. 

We note that due to the lack of quantitative theories, the fitting formula used here is 
somewhat empirical, and certainly could be replaced by more sophisticated formulae with 
different number of components. However, the essence of SWT or any other equivalent is to 
capture the anomalous rapid spectral weight growth below at the gap energy scale beyond 
the Fermi-liquid-based descriptions and to offer a direct measure of the superconducting or- 
der parameter which, in return, facilitates examination of various theoretical models. There 
have been a few theoretical proposals of such a quasiparticle peak created in the supercon- 
ducting state by spectral weight transfer [2^, and its weight was related to the condensate 
fraction or ODLRO. However, most of these theories are based on underdoped cuprates, 
while the overdoped regime is seldom addressed. The observed correlation between the 
SWT and superfiuid density throughout the entire phase diagram calls for further theoreti- 
cal efforts. 
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Energy relative to Ep (eV) 

FIG. 1: ARPES spectra for optimally doped Bi2212 as a function of a) T and b) momentum in 
the superconducting state along (0.457r, 0) — (7r,0). Features in the circled region are due to the 
superstructures, c) The extracted dispersions of the peak and hump based on the spectral maxima 
of features, d) The calculated peak-hump spacing is compared with the theoretical prediction for 
bilayer energy splitting in the superconducting state. 
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Energy relative to Ep (meV) 

FIG. 2: T-dependent ARPES spectra at (vr, 0) of OD65 in (a) fitted by 2-component model, (b) 
by 3-component model including /iz^-dependent spectra on Pb-OD65 at (c) T=10 K and (d) T=80 
K (see text). The AB hump lies lower in energy than the BB hump. The width difference of the 
SCP is due to resolution variations from 10 to 18 meV at different hi/^s. 
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FIG. 3: The fitted a) T-dependent energy positions of the maxima of various features and (b) 
/iz^-dependent matrix elements of different spectral features. 




FIG. 4: (a) The T-dependent SWT's given by 3-component and 4-component fit are shown. The 
interpolations are based on a trial function SWT = SWTq[1 - (r/Tc)^]Q on the T-dependence 
of SWT (averaged) . The fitted value, C, ~ 2.5, agrees well with the T-dependent /iSR result 
on Tl2Ba2CuO(i+s at the same doping level (OD66 T12201) (after Ref.[18]). (b) The doping- 
dependence of SPR at hv=22.7 eV supplemented by the global fit result on OD65. 
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